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We perform a lattice mass analysis in the flavour sin- 
glet pseudoscalar channel on the SESAM and T%L full QCD 
vacuum configurations, with 2 active flavours of dynamical 
Wilson fermions at (3 = 5.6. At our inverse lattice spac- 
ing, a -1 ~ 2.3 GeV, we retrieve by a chiral extrapolation to 
the physical light quark masses the value my = S.T^m^. 
A crude extrapolation from (Nf = 3) phenomenology would 
suggest m n i ~ 5Am 7r for Nf = 2 QCD. We verify that the 
mass gap between the singlet state n' and the n flavour triplet 
state is due to gauge configurations with non-trivial topology. 



I. INTRODUCTION 

Lattice gauge theory (LGT) has been established as 
the standard method to deal with infrared aspects of 
quantum chromodynamics (QCD). Recently, the light 
hadronic flavour non-singlet masses have been accurately 
determined by Nf = 2 QCD simulations on the teracom- 
puting scale jjj . Unfortunately, the situation is much less 
clear when it comes to the interesting physics of flavour 
symmetric hadronic states (like the ?/-meson) which are 
expected to be influenced by the topological properties of 
the QCD vacuum: we have to await multi-teracomputing 
to see it settled. 

The problem is due to the very occurrence of Zweig- 
rule forbidden contributions to the singlet hadronic prop- 
agator in form of disconnected diagrams, as already 
known from early feasibility studies in quenched QCD M. 
The difficulties arise for three reasons: (a) eo ipso, such 
disconnected correlators induce a high level of gauge field 
noise into the calculations; (b) their computation is costly 
as it involves momentum zero projections of quark loops, 
the evaluation of which requires the use of stochastic esti- 
mator techniques; (c) the propagator of a flavour singlet 
pseudoscalar meson, C n i , turns out to be the difference 
between connected and disconnected diagrams with the 
possibility of numerical cancellations. 

Indeed, at large Euclidean time separations, t, these 
cancellations are doomed to be strong if they are to ren- 



der the large empirical flavour singlet/non-singlet mass 
gap [], Mq = Mi — M%, As a consequence, however, the 
signal-to-noise ratio becomes a serious problem for the di- 
rect lattice approach to flavour singlet objects and makes 
it hard to keep control on systematic errors. For all these 
reasons, ab initio full QCD lattice investigations of the n' 
mass have not yet overcome an exploratory stage. There 
is of course a way to avoid all this by following an indirect 
strategy and taking resort to the assumptions underly- 
ing the Witten-Veneziano formula ||; this workaround 
amounts to determining the mass gap from quenched lat- 
tice determinations of the topological susceptibility 0]. 

Previous pioneering work in full QCD largely focussed 
on a two-step recipe to deal with the above problems: 
(i) determine at large t and (ii) compute the mass 
gap 0, mo, from the ratio of connected and disconnected 
correlators, R(t) = C dl sc{t)/C conn (t) (| [J in the range of 
smallish t- values. In our present approach we seek for a t- 
window within which a straightforward (one-step) flavour 
singlet propagator analysis can be pertinently achieved. 

In a recent study we already applied improved stochas- 
tic estimator techniques — as geared previously for cop- 
ing with disconnected operator insertions in the context 
of hadronic matrix elements || — on the flavour singlet 
correlators, with pointlike sources §JlO). In this paper 
we shall show by a mass plateau analysis that a standard 
mass computation on the flavour singlet propagator itself 
will become feasible with reasonable control of systematic 
errors, once smeared operators are used. 



1 ln our Nf — 2 world we have a triplet (rather than an 
octet) of flavour non-singlet mesons. Moreover, working with 
mass-degenerate quarks, our 7r's are exactly mass degenerate 
too. 

2 Upper (lower) case letters refer to masses in physical (lat- 
tice) units. 
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II. LATTICE PREREQUISITES 

We consider the pseudoscalar flavour singlet operator 
in a flavour symmetric theory 



S(x) = ^g i (x)7 5 %(x) 

i=l 



(1) 



with Nf flavours. By the usual Wick contraction it leads 
to the flavour singlet propagator in terms of the inverse 
Dirac operator, A = D~ x : 



C n >(0\x) ~ (Nftr(A(0\x)A^(Q\x)) 

-N 2 f tv( l5 AH0\0))tv( l5 A(x\x))) , 



(2) 



which is a sum of fermionic connected and disconnected 
contributions with traces to be taken in the spin and 
colour spaces. In the rest of the paper we shall refer to 
them as 'one-loop' and 'two-loop' contributions, respec- 
tively. The traces are computed with Z2 noise sources 
including diagonal improvement as explained in ref. |9| . 
The momentum zero projection 



C v ,(t) = (S(t)S(0)) 

c 



(S(t)S(0)) di . 



(3) 



is expected to decay exponentially, ~ exp(— m v >t), and 
thus to reveal the flavour singlet mass, m, ; < . On a toroidal 
lattice with temporal extent T one should encounter the 
usual cosh behaviour at large values of t and T — t 



C v >(t) — ► exp(— m v 't) + exp(— m v >(T — t)) 



(4) 



From this parametrization effective masses, mL , can be 
retrieved by solving the implicit equations 

CV(t+l) _ e X p(-m*,(i + l))+e X p(-m*,(T-t-l)) 



exp(— m*,i) + exp(-m*,(T - t)) 



(5) 



For sufficiently large values of t, the effective masses 
should saturate into a plateau. The crucial question is, 
however, whether one can establish a t-window of obser- 
vation that reveals a definite plateau behaviour of m*, 
before noise takes over. 



A. Operator smearing 

The building blocks for hadronic observables are the 
quark propagators, £, which may be computed by solving 
the (discretized) Dirac equation with appropriate source 
vectors, <f> s (z), on the lattice: 



D(z,x)£(x) = 4> s (z) 



(6) 



In standard spectrum analysis one generally applies some 
kind of spatial smearing to the hadron source (located at 
t = 0) in order to enhance the ground state signals of 



TABLE I. Simulation parameters used at /? = 5.6 and 
numbers of stochastic sources (used with local and smeared 
operators, N<? st , N^). Last column: numbers of available 
decorrelated vacuum field configurations, Nconf- 
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195 


0.1575 
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*32 


400 


400 


195 


0.1575 
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400 
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156 


0.1580 
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24 :l 


*40 


100 


100 


156 



the resulting hadronic operators at medium values of t. 
Needless to say, this appears to be all the more necessary 
in the present context where - as explained above - we 
are faced both with (a) cancellations (between C conn and 
Cdisc) and (b) noisier signals (from C disc ). 

We used our smearing procedure as applied in the anal- 
ysis of light non-singlet masses |ll[] ; it is characterized by 
N diffusive iteration steps 



>(*) = 



1 



l + 6cv 



(x,t) +a^2^(x + 



(7) 



where the index 'p.t.' stands for 'parallel transported', 
and the sum extends over the six spatial neighbours of 
x. For c^i ' we start out from pointlike sources for the 
connected and ^2-noise nonlocal sources for the discon- 
nected diagrams. In this way the bilinear quark opera- 
tors, Eq. (|l|), were computed after N = 25 such smearing 
steps, with the value a — 4.0. The smearing procedure 
was applied to meson sources as well as to sinks - both 
for Cdisc and C conn , in order to correctly maintain their 
relative normalisations. 

In table | we list the run parameters of our simula- 
tions, which make use of vacuum field configurations gen- 
erated by the SESAM (16 3 x 32 lattice @) and the T^L 
(24 3 x 40 lattice |]|) collaborations, both with Nf = 2 
and /3 = 5.6. We have used five different sea quark masses 
and two different lattice sizes to gain some control on 
finite-size effects. While the number of vacuum configu- 
rations varies from 156 to 195, the number of independent 
stochastic sources has been chosen to be 400 on the small 
lattices, both for local (It) and smeared (srri) operators. 
On the large lattices 100 (400 for n sea = 0.1575-ZZ) source 
vectors were used. 

Fig. illustrates the quality of our data in terms of 
the one-loop and two-loop correlators at the lightest sea 
quark mass on the SESAM lattice, computed with point- 
like (upper figure) and smeared operators (lower figure), 
using N est = 400 stochastic Z2-noise sources with diag- 
onal improvement Q). The errors quoted are statistical 
and have been obtained by jackknifmg. We find a marked 
improvement of the signal by help of source smearing in 
the regime 5 < t < 12. 
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FIG. 1. Effect of smearing on the correlation functions. 
SESAM configurations at the lightest sea quark mass, 
Ksea — 0.1575. Top figure: local operator; bottom figure: 
smeared operator. Upper (lower) data set refers to one-loop 
(two-loop) contributions. 

One should remember that the data points in Fig. |l| 
suffer from two kinds of stochastic noise: the one from 
the gauge fields and the one from the noisy sources. In 
order to disentangle them it is highly instructive to study 
the error on the two-loop signal, a, as a function of the 
number of stochastic sources, N est . In Fig. |2| we plot 
this quantity for the SESAM ensemble at n sea = 0.1570 
on a time slice of interest, t = 8. At sufficiently large 
values of N est the parametrization 



-'conf 



(8) 



is expected to describe the superposition of errors from 
source and gauge field fluctuations. Therefore, the er- 
ror analysis can provide a useful check on the quality of 
the stochastic estimator outputs. The data in Fig. g in- 
deed yield convincing evidence for early asymptotic N est - 
dependence, with a threshold value N est ~ 64. Moreover, 
we find that on the SESAM sample the genuine gauge 
field noise (as indicated by the horizontal asymptotic 




iV es t — * oo line) prevails once we choose N est > 100. As 
to the subasymptotic regime, one might attribute the ap- 
parent 'non-standard' behaviour of a to pollutions from 
subleading, non-trace terms in the stochastic estimate of 
the loops. 

The main message from Fig. ^| is that the Z2 noise 
method does provide reasonable parametric control over 
the additional fluctuations induced by the stochastic es- 
timator on the observable. This control cannot be taken 
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FIG. 2. Error, a, of two-loop signals versus N eat at 
Ksea = 0.1570 on time slice t = 8, for smeared sources and 
sinks. The curve is the best fit according to the parametrisa- 
tion given by Eq. (^). 

for granted when applying the volume source technique 
(12 color-spin explicit inversions per configuration) which 
refrains from using stochastic sources and relies fully on 
gauge invariance and gauge noise for the suppression of 
nondiagonal contributions to the trace estimates . For 
comparison, however, we have included the correspond- 
ing error as we computed it on our gauge field ensemble. 
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FIG. 3. t-dependence of the two-loop signal, on a single 
configuration with smeared operator, for various values of 
N est . 

Complementary to these considerations one may study 
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FIG. 4. Ground state dominance of the r/-propagator 
i^sea = 0.1580), with smeared sources and sinks. Top: with 
local operators, bottom: with smeared operators. 



the overall (in t) effects of finite source sampling on the 
estimate of the two-loop correlator for a particular gauge 
configuration. Fig. |] illustrates, again at n sea = 0.1570 
on the small lattice, the kind of fluctuations induced by 
the stochastic sources of the two-loop correlator with 
smeared operators at various values of N est which we ran 
up to 2048 in this case. It appears that on our sample 
the Zi noise injected from the sources into the correla- 
tor is adequately suppressed at N est « 400. This justifies 
again that on the small lattices N est — 400 is a reasonable 
choice for the present study. On the large lattices, how- 
ever, enhanced self-averaging effects allow for a smaller 
number of stochastic sources, N est ~ 100. 

We are now in the position to study plateau formation 
on our ensemble of vacuum configurations. 



B. Plateaus of effective masses from smearing 

The effect of smearing on the flavour symmetric cor- 
relator, Crji , is visualized in the comparative twin plot 
of Fig. ||, as obtained at our smallest quark mass on the 
T%L lattice. Prima vista we do find reasonable signals 



on this correlator up to t w 10. Moreover, by inspection 
of the cosh-fits, we find a considerable decrease of excited 
state contributions as a result of smearing. 
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FIG. 5. Effective masses in the flavour singlet and 'octet' 
channels, with pointlike (top) and smeared operators (bot- 
tom), at K sea = 0.1570. 

Let us now scrutinize the situation by turning to the 
analysis of effective masses as extracted from Eq. (||) . Wc 
have seen in ref. that, given our sample sizes, the use 
of local sources and sinks does not provide sufficient res- 
olution to reveal such plateau formation in the effective 
flavour singlet mass plots. This situation is illustrated 
in Fig. ^| where we confront, at a particular intermediate 



TABLE II. Fit ranges 
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FIG. 6. Plateau formation in the effective n'- and n masses 
with smeared operators at various sea quark masses. At 
Ksea = 0.1575, ^'-results from SESAM and T^L configura- 
tions are plotted separately. 
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FIG. 7. Chiral extrapolations of m v , (top figure) and m,/ 
linear in quark mass. 



the singlet channel. 

In Fig. U we display the evidence for plateau formation 
through operator smearing for the remaining sea quark 
masses in the range 0.1565 < K sea < 0.1580. Here again 
the 'octet' channel masses are included for reference in 
order to enable judgement on the sensitivity for mass gap 
determinations. We emphasize again that all singlet data 
are obtained after symmetric source and sink smearing, 
as described above. It appears that smearing meets the 
expectation by increasing the ground state overlap: this 
opens the window of observation for a mass plateau from 
t = 5 onwards where statistical errors are still 'tolerable'. 
By comparing the SESAM and T^L data sets at n sea — 
0.1575 we find no evidence for a volume effect on 



sea quark mass (K sea — 0.1570), effective pseudoscalar 
masses obtained both with and without operator smear- 
ing. Obviously, with pointlike operators, one has to re- 
sort to bona fide single cosh fits on the correlators with- 
out any kind of systematic error control on the extracted 
flavour singlet 'masses'. After source and sink smearing, 
however, our data begins to reveal plateau formation in 



III. PHYSICS ANALYSIS 

Encouraged by the apparent plateau formation we pro- 
ceed next to carry out mass fits based on a single cosh 
ansatz, with rj-ranges as listed in table |J. F° r reference 
we have also included information about the fit-ranges 
previously used with local sources M. 
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A. Chiral extrapolations 

Because of the well-known technical limitations of the 
hybrid Monte Carlo algorithm JTs| the SESAM and TxL 
configurations correspond to two mass-degenerate light 
sea quark flavours (Nf — 2), with the unrenormalizcd 
mass value 



m q = M q a = 1/2(k 1 - k" 1 ) . 



(9) 



From our previous light spectrum analysis |Tl| ] we quote 
the lattice spacing 

a; 1 (^ 3 , lt ) = 2.302(64)GeV (10) 

and the critical and physical light quark n values: 

k c = 0.158507(44) , K Ught = 0.158462(42) . (11) 

Our data do not allow to decide whether it is m 2 , or m v > 
that follows a linear quark mass dependence: as shown 
in Fig. 0, both ansatze render x 2 /d-o.f. ~ O(l). We 
emphasize in this context that we make no distinction 
between sea and valence quarks as we choose the quark 
masses in the fermion loops to equal the sea quark masses 
(symmetric extrapolation in the sense of ref. |ll| ) . 

We display the results on the rf mass and the mass 
gap 

ml = ml, - ml , (12) 



for both forms of extrapolation in table III, For com- 
parison, we have also included previous estimates as ob- 
tained by using local sources 



B. Comparison to experiment 

In the Nf = 2 world of our simulations, according to 
Eq. (Q) , we would not expect to encounter the full effect of 
Zweig rule forbidden diagrams, and hence we anticipate 
to underestimate the real world rj mass, m oxp , (plotted 
as open squares in Fig. 0). 

From the experimental mass splitting 



Mo,n F =3 



lvl ri',n F =3 



M 



8 > 



M 2 = 2M 2 K -M 2 (13) 



we therefore compute, in the spirit of the Witten- 
Veneziano formula, 



TABLE III. n' and m results 
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m n i 
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.245(40) 


551(85) 


565(92) 


Nf 
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FIG. 8. Chiral extrapolation of mo. 



Ml = 2N fXq /K . (14) 
the 'pseudoexperimental' value, M 8 , in our Nf = 2 world: 

M 2 S = 2/3 Ml,, , + Ml = (715 MeV) 2 . (15) 

This value corresponds in lattice units to the full squares 
marked 'm s ' in the two alternative chiral extrapolations 
shown in Fig. |t]. Let us compare this latter value with 
the lattice Nf = 2-prediction in terms of numbers: when 
we set the scale by the p mass, the extrapolation of m 2 , 
to the physical quark mass yields the value 



M 2 , 

rj' 



(520^ 25 MeV) 2 



(16) 



at our lattice spacing. We have linearly added the differ- 
ence between the two extrapolations (of m v > and m 2 , ) to 
the statistical error to accomodate systematic uncertain- 
ties. The difference between this value and the 715 MeV 
of Eq. JTe| ) may be attributed to finite-a effects. For 
comparison, we also performed a linear extrapolation of 
the mass splitting, mo, as given in Eq. ( |l2] ) which shows 
only little quark mass dependency (see Fig. @). This is 
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FIG. 9. Ratio of Cdiac/C CO nn for K sea = .1575 (top figure) 
with cuts in topological charge as explained in the text. 
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FIG. 10. Effective rj (top) and ■k (bottom) masses for 
K sea = .1575, with cuts in topological charge as explained in 
the text. 



consistent with the weak dependence of F v on the quark 
mass observed on SES AM configurations (jllj . From the 
value F n = 116(8) MeV obtained for N F = 2 QCD at 
our lattice spacing and Eq. ( |l6| ) we obtain the estimate 
X q = (170^^5 MeVY for the quenched topological sus- 
ceptibility of Eq. (|14|) , which is consistent with the value 
180 MeV from 3 flavour phenomenology |B[. 



C. Impact of topology 

The Witten-Veneziano mass formula, Eq. (jlj), relates 
the difference between the r]' mass and the flavour non- 
singlet pseudoscalar mass to the topological susceptibil- 
ity Xq °f t ne quenched gauge vacuum. This scenario 
motivates us to investigate in full QCD whether the ra- 
tio R,Q(t) = Cdi SC (t)Q/C conn (t)Q, whose deviations from 
zero give rise to the observed mass gap, is correlated with 
|Q|, the modulus of the topological charge, configura- 
tion by configuration. If wc restricted our analysis for 
instance to gauge configurations with Q = only, the 
topological susceptibility, x = (Q 2 )/V, determined on 
this subensemble, would vanish as well and we might ex- 



pect the 7T to be mass degenerate with the rj . On the 
other hand, if we rejected configurations with small \Q\- 
values, the effective x ° n ths remaining sample would be 
enhanced and the generated mass gap increased. 

In Fig. |^ we show, for n sea = 0.1575 and the small 
lattice, the quantity Rc){t) with cuts applied according 
to \Q\ < 1.5 (topOl) and \Q\ > 1.5 (top24), the topologi- 
cal charge being determined as in ref. |16|. The value of 
1.5 was chosen such as to obtain two ensembles of com- 
parable statistics in topologically different vacuum sec- 
tors. We do find a definite dependence of Rq on |Q|. 
Note in particular that the disconnected piece vanishes 
in the vacuum sector with small values of |Q|! This fea- 
ture reflects itself of course in the corresponding effec- 
tive masses of the flavour singlet and non-singlet mesons. 
This is shown in Fig. HG: the restricted flavour singlet 



v' I|Q|<i.5j 



turns out to be identical to the 'octet' 



mass, m 

meson mass — m w . Accordingly, the flavour singlet- 
non-singlet mass gap is due to nontrivial topological vac- 
uum structures. 

On the other hand, the 'octet' meson mass appears to 
be not at all sensitive to such restrictions to topological 
sectors; this seems to be a general feature of flavour non- 
singlet light hadron spectrum observables Il7| . 



IV. SUMMARY AND CONCLUSIONS 

Using smeared operators and reasonable source statis- 
tics on the SESAM and T^L samples of QCD vacuum 
configurations we found clear indications of plateau for- 
mation in the effective flavour singlet pseudoscalar mass 
plot in the intermediate i-regime. The rj mass is defi- 
nitely sensitive to the topological structure of the QCD 
vacuum. In our two-flavour simulation its actual value 
after chiral extrapolation turns out to be in qualitative 
agreement with the expectation from experiment, but 
further studies are needed to pin down finite-a effects. 

At this stage the statistical errors on the singlet masses 
are mostly due to gauge field fluctations and by a factor 
~ 5 larger than for the non-singlet ones. Clearly, the 
next generation teracomputers will open the door to go 
for lattice determinations of Zweig-rulc forbidden objects 
with an accuracy known so far only from light non-singlct 
hadron spectroscopy. 

In the meantime we are working on computational 
techniques to determine quark loops in the regime of 
quark masses lighter than attained so far C3] . 
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